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Abstract 

We derive different representations of compact QED fixed to Lan- 
dau gauge by the lattice Faddeev-Popov procedure. Our analysis finds 
that (A)Nielsen-Olesen vortices arising from the compactness of the 
gauge-fixing action are quenched, that is, the Faddeev-Popov determi- 
nant cancels them out and they do not influence correlation functions 
such as the photon propagator; (B)Dirac strings are responsible for 
the nonzero mass pole of the photon propagator. Since in D = 3 + 1 
the photon mass undergoes a rapid drop to zero at /3 C , the deconfine- 
ment point, this result predicts that Dirac strings must be sufficiently 
dilute at (3 > (3 C . Indeed, numerical simulations reveal that the string 
density undergoes a rapid drop to near zero at (3 ~ (3 C . 



1 Introduction and Results 



Gauge fixing is essential to several potentially physically relevant lattice com- 
putations, of which we mention two: (i)Partial gauge fixing of SU (3) to resid- 
ual U{Y) X U{Y) is necessary to define abelian projection monopoles [1, 2] 
whose currents, as a working hypothesis, may be the underlying confinement 
mechanism of QCD. (ii)Computing fixed-gauge lattice matrix elements may 
be a way to determine continuum-lattice renormalization or "matching" coef- 
ficients. In particular, such coefficients are necessary (in a certain approach) 
for trying to exhibit the Delta 1=1/2 Rule on the lattice [3]. 

With such ultimate motivations, there have been many numerical stud- 
ies of lattice gauge fixing and the gauge dependence of such gauge variant 
quantities as the abelian projection monopole density [1, 2, 4] and effective 
gluon, quark, photon and electron masses [5, 6, 7]. In addition, the gauge de- 
pendence of quark masses has been analytically computed for certain gauges 
in the strong-coupling expansion [8] and in the Schwinger model [9]. 

While compact QED("cQED") in the strong coupling regime in the ab- 
sence of gauge fixing is a well-studied model of a confining gauge theory [11], 
its fixed-gauge features are less known and also nontrivial [6, 12]. In this pa- 
per, we report on a numerical and analytical study of cQED fixed to Landau 
gauge. As described in Section 2 both the photon mass 1 and Dirac string den- 
sity("kink" density) drop dramatically from nonzero to near zero at f3 ~ /3 C , 
the deconfinement point in D = 3 + 1 dimensions. Is this a coincidence or are 
Dirac strings dynamically related to photon mass? Sections 3 and 4 cast the 

1 We stress that photon "mass" in this paper refers to the pole of the photon propagator. 
Photon mass thusly defined is gauge variant and is not obviously related to such physical 
length scales as the electric penetration depth, which is gauge invariant. 
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lattice Faddeev-Popov procedure into lattice differential forms notation [13], 
using which we show that: 

• (A)Abrikosov-Nielsen-Olesen [15] vortices from the lattice gauge fixing 
"spin glass" action are quenched and, so, do not contribute to photon 
mass; 

• (B)Disorder caused by Dirac strings is responsible for photon mass [12]. 

2 Numerical Results 

The cQED action is 

S c = /3 "£(1- cos F^) (1) 

/i<i/ 

where i^fA] = d ll A v — d^A^ and £ [— 7r, 7r). While S c is gauge invariant, 

F^ = 0^ + 2* (2) 

decomposes into a gauge invariant electromagnetic field part £ [— 7r,7r) 
and an integral part N^. If / on a plaquette, we say that the pla- 
quette has a "kink." Dirac strings, which costs zero action, and magnetic 
monopoles, the Dirac string endpoints which cost action, are comprised of 
kinks. Since the condition A^ £ [— 7r,7r) is enforced by 2tt modding, 
transforms nontrivially under local gauge transformations which push A^ 
outside of [— 7r,7r) and, so, Dirac strings are gauge variant. Monopoles are 
gauge invariant. 

Landau gauge in this paper is defined as in Ref. [6]. As described in 
Appendix A, an alternative definition — henceforth called "cLandau" gauge — 
was adopted in Ref. [12]. cLandau gauge is a compactness-preserving gauge 
which can heuristically be thought of as (9 |U A |U )mod(27r) = 0. In the ter- 
minology of this paper, the "Landau" gauge photon mass values reported in 
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[12] are really in cLandau (not Landau) gauge. Photon masses quoted in this 
paper are in Landau gauge. 

While the gauge variant Dirac strings do not have any direct physical 
effects, they play a central role in the fixed-gauge sector. Figure 1 depicts a 
typical Landau gauge configuration in cQED 2 +i on a 25 3 lattice with periodic 
boundary conditions. The thin lines mark the Dirac strings, the big dots the 
monopoles and antimonopoles. The reader should be convinced that Dirac 
strings either connect monopole-antimonopoles pairs or form closed loops. 
In Landau gauge the links tend to U x ^ — > 1 at large f3 and at f3 = 2.3 the 
operator 

1 3 

F ( x ) = 9 J2\- cosA v( x ) + cosA ii( x - A)] (3) 

has a typical value of F(x) ~ 0.95. However, 3% of the sites obey F(x) < 0.77. 
These "small-F" sites are indicated in Figure 1 by the small dots. 

Let us mention that we have seen analogous small-F sites in the course 
of gauge fixing SU (3) gauge theory to maximal abelian gauge and Landau 
gauge. While we do not claim to understand at this point that small-F sites of 
cQED have anything in common with QCD small-F sites, we have determined 
that SU{3) small-F sites have (quite dramatically) fractal dimension Dj ~ 2 
for Landau gauge and, less dramatically, a bit smaller Dj for maximal abelian 
gauge. These SU (3) results are preliminary and will not be further discussed 
here [2]. 

In Figure 1 the small-F sites cluster around closed and open Dirac strings 
because the vector potential around a Dirac string in the z direction is (in con- 
tinuum cylindrical coordinates) A strmg = (j> , 1 and around a monopole 

V x 2 +y 2 

at x = is (in spherical coordinates) A mono = <f> 1 ~ ;° s / • Consequently, near 
Dirac strings and monopoles the smoothness of the vector potential is dis- 
rupted. This disruption disorders the fixed-gauge sector of the cQED vac- 
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Figure 1: A typical f3 = 2.3 cQED 2 +i gauge configuration on a 25 3 periodic 
lattice after about ~ 2000 30-hit Metropolis/pseudo-heatbath thermalization 
sweeps. The configuration is fixed to Landau gauge. The thin lines mark the 
Dirac strings, the bigger dots the magnetic monopoles and antimonopoles, 
and the smaller dots the locations of the smallest 3% of F(x). 
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uum. As described in Section 4, one consequence of this disorder is that the 
mass pole M 1 of the photon propagator is nontrivial. 

Figure 2 depicts the Landau gauge kink number density 

_ ^kinks in the lattice 
^plaquettes in the lattice 

as a function of f3 in D = 3 + 1 dimensions. Since is gauge variant, p^ 
varies with gauge. In Landau gauge, like the photon mass reported in Ref. [6] 
Pn drops dramatically to near zero at f3 ~ f3 c . Figures 3 and 4 show that 
photon mass (when it is nonzero) is, approximately, linearly correlated with 
Pn in both D = 3 + 1 and D = 2 + 1 dimensions. Note that in D = 2 + 1 
we are quoting the dimensionless quantity @M 1 rather than M 1 . The lines 
in Figures 3 and 4 are to guide-the-eye. 

In Landau gauge p^ and monopole density pu are positively correlated. 
(pM also drops sharply near f3 c [16].) This is because Landau gauge suppresses 
Dirac string length — most of the contribution to p^ comes from unavoidable 
Dirac strings which connect monopoles-antimonopole pairs. Moreover, while 
monopole-prohibition (the procedure described in Ref. [12]) decreases both 
Pn and M 1 substantially in Landau gauge, p^ and M 1 decrease together so 
that they remain on the original p^-M^ line (see Figure 4). 

On the other hand, since (9 |U A |U )mod(27r) = leaves ±7r ambiguities on 

cLandau gauge does not suppress Dirac strings. In this gauge monopole 
prohibition does not decrease kink density much — most of the kinks are in 
closed Dirac strings. Hence in cLandau gauge monopole prohibition does not 
significantly change the photon mass [12]. This is another piece of evidence 
linking photon mass to Dirac string density. 

In summary, Landau gauge suppresses Dirac strings whereas other gauges 
such as cLandau do not. In any gauge and in both D = 2 + 1 and D = 3 + 1, 
when Dirac strings are dilute the photon mass is small; as Dirac strings be- 
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Figure 2: The Landau gauge kink number density p^ as a function of f3 in 
cQED 3+ i on an 8 3 X 16 periodic lattice. The first configurations at each 
f3 is thermalized starting from cold-start by 2000 30-hit Metropolis/pseudo- 
heatbath sweeps. The error bars indicate the standard deviation of p^ in 20 
successive gauge configurations separated by 40 sweeps. 
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Figure 3: A scatter plot showing the correlation between the Landau gauge 
photon mass M 1 and the Landau gauge kink number density pjy at common 
/3s in cQED 3+ i on an 8 3 X 16 lattice. The gauge configuration generation 
is done as for Figure 2 except that in this case 500 gauge configurations 
separated by 5 sweeps were used for each data point. We stress that each 
data point in this Figure is generated fresh from a cold-start independently 
of any other data point so there is statistical correlation between the plotted 
points. Errors are jackknife errors. The line is to guide-the-eye. 
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Figure 4: Same as Figure 3 except (i)we plot f3 X M 1 rather than Af 7 ; 
(ii)the calculation is in cQED 2 +i rather than in cQED 3+ i; (iii)we include 
data points(the asterisks) for which monopoles are prohibited during the 
Metropolis/pseudo-heatbath updating sweeps. 




come denser, the photon mass becomes larger. This monotonic relationship 
holds whether monopoles are prohibited or not. 

3 Quenching of Spin Glass Vortices 

In this Section we apply lattice differential form and BKT [10] transforma- 
tion ideas, previously exploited to expose cQED without gauge fixing [11] 
and the compact abelian Higgs model [13], to analyze the lattice Faddeev- 
Popov determinant arising in lattice gauge fixing. The strategy is to perform 
a variables change so that the underlying excitations — spin glass vortices 
arising from the periodicity of the lattice gauge fixing action — are explicitly 
displayed. We will show that these vortices are quenched by the Faddeev- 
Popov determinant and, so, do not play an important role in disordering 
such gauge variant correlation functions as the photon propagator. These 
spin glass vortices should not be confused with Dirac strings, which come 
from the cQED action. 

The gauge fixed expectation value of a gauge variant operator O(A) ^ 
0(A — d(p) as computed in numerical simulations is 

JVA exp{-S c (A)} 0(A axed) 

gauge / K \ 

JVA' exp{-S c (A')} { } 

where path integral measure DA is not gauge-fixed. Definition (5) is equiv- 
alent to the Faddeev-Popov approach since, with a suitable choice of gauge 
fixing action S'gf(A) ^ S g {(A — d(p) } operator 0(Aaxed ) can be represented 

gauge 

as 

0(A axed) = A FP (A) fvipexp{-S g {{A-dip)}0{A-dip). (6) 

gauge J 

The substitution (9 = 1 into (6) yields for the inverse Faddeev-Popov deter- 
minant 

Ap^A) = f XV exp {-S^A -dtp)}. (7) 

J — IT 



CM fixed 
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Note that App(A) = App{A — dip). Hence as in continuum models the 
Faddeev-Popov determinant is gauge invariant. Inserting (6) into (5) and us- 
ing the gauge invariance ol S c and App to transform away the (^-dependence 
ol the DA integrand yields [7, 8] 

(0(A^))=Z[0}/Z[1] , (8) 

gauge 

Z[0] = JVA exp{-S c (A)} G(A) 0(A), (9) 

G(A) = App(A) exp {S g{ (A)} . (10) 

While our numerical simulations are performed in the Wilson represen- 
tation, we will now adopt the Villain representation, which is more amenable 
to formal analysis. There is no doubt that for our purposes the Wilson and 
Villain approaches are qualitatively (albeit not quantitatively) equivalent. In 
the Villain representation 

5 c (A) = -log J2 exp{-|||dA-27riV|| 2 }, (11) 
jv(c 2 )g^ I z ) 

S si (A) = -log £ exp{-^||A- 27 r/|| 2 }. (12) 

The k — > oo limit corresponds to minimizing S L of Eq. (A.l). Following from 
(7) and (12), the inverse Faddeev-Popov determinant is 

^fp(A)= £ flV exp{-^||A-d^-27r/|| 2 }. (13) 

In (13) App is equivalent to the partition function of the compact abelian 
Higgs model in an external electromagnetic field. 

One can manipulate lattice differential forms just like continuum ones. 
Our integration by parts convention is (<p,5ip) = -\-(d<p,ip) so that 8 = + * 
d* in any dimensions [13]. The Laplacian is A = Sd + dS and j-forms are 
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generically denoted c r Zero modes, which are related to Gribov copies, will 
be ignored in this paper. The Hodge-DeRham decomposition of any real- 
valued form A follows from 

1 = 8A~M + dA _1 6 (14) 

which implies A = 6A~ 1 F + dco where dF = = 8u. For integer j-form /, 
the Hodge-DeRham decomposition takes the slightly modified form [13] 

l(cj) = 5A _1 u(c j+ i) + dA- 1 r(c J _ 1 ), v = dl, dv = 0, 8r = 0. (15) 

Then one can show [14] that up to an overall volume proportionality factor 

E <x E E (i6) 

l(cj)E!Z {v(c j+1 )E!Z\dv=0} {r(c J _ 1 )e^|5r=0} 

Hodge-DeRham decomposition of A and /, application of (16), and ab- 
sorption of all the exact pieces into dtp yields 

Vif exp |- - (dip, dip) + K (<p, 8A)j (17) 

where A(A) = exp | — 1|| A|| 2 | and 

Z V (A) = E exp ^—2tt 2 k(v, A~ 1 v) + 2ttk(v, A _1 dA)| . (18) 

{v(c 2 )£Z\dv=0} 

In the RHS of (17) the overall proportionality constant arising from (16), 
which has no dynamical consequences, has been dropped. 

Due to the condition 8*v = 0(dv = 0), the vortices *v in Z v are closed 
Abrikosov strings in D = 2 + 1 and closed Nielsen-Olesen string worldsheets 
in D = 3 + 1. Due to the factorization of A and Z v from the p integral, 
electromagnetic sources from A and vortices from Z v are quenched in Z[0], 
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that is, A and Z v drop out of G(A) in Eqs. (9-10). Putting all the pieces 
together yields 

G(A) = E r( c )e^exp {-27r 2 K(r, A~V) + 2 7 r«(A- 1 r, 8 A)} 

Ho *V ex P {-f( d ^ <V) + <^)} 

S exp {-2^4(r - M), A^(r - M))) + f (M, A~HA)} 

_ r{c )E% 

exp{f(<5A, A^SA)} 
= £ ^{-2* 2 "((r-±6A),A-\r-±6A))} (19) 

where 8r(c ) = is automatic since r is a 0-form. Taking the n — > 00 limit 
yields the lattice Landau gauge gauge fixing condition 

lira G(A)= ]T 6(SA-2wr). (20) 

Eq. (20) implies that 8 A is ambiguous up to 2irr where r(c ) G 8r = 0. 
These integer gauge ambiguities do not affect Dirac strings or any operators 
which are periodic in A^. 

4 Origin of Photon Mass 

The photon mass contribution from the gauge fixing action and the Faddeev- 
Popov determinant is given by the ||A|| 2 term in the functional Taylor ex- 
pansion of log G(A). Since log G depends only on the longitudinal part of A, 
it cannot have an ||A|| 2 term in its Taylor expansion and, hence, gauge fixing 
operator G(A) cannot generate a nonzero photon propagator mass pole. The 
photon mass must come from S c . 

To isolate where the nonzero M 1 comes from, introduce the noncompact 
QED action 

S nc (A) = \ ||dA|| 2 (21) 
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where A £ ( — 00, 00). S nc is just the lattice transcription ol massless contin- 
uum QED whose Landau gauge photon propagator is massless. Expression 
(8) can be rewritten as 

Z[0] = J VA F(A) exp- 5 " c(A) G(A) 0(A), (22) 

F(A) = exp~ s ^ exp +5 " c ^ . (23) 

Since S nc and log G have no photon mass contribution, the photon propagator 
pole is determined completely by 

\ogF{A) = - l -M 2 1 || A|| 2 + ••• (24) 

To probe into the origin of Af 7 , make the variables change 

N = <5A _1 m + dA" 1 ^, dm = 0, 8q = (25) 

on N of Eq. (11). (D — 3)-form *m is the monopole current and, since 

S*N = *m, (26) 

*N is the Dirac "sheet" in D = 3 + 1 dimensions. Application of Eq. (16) 
yields 

F(A) = C m exp{-27r 2 /3( g ,A- 1 g ) + 27r/3(A, (Z )} (27) 

{g(ci)G^|5g=0} 

where 

C m = exp{-27r 2 /3(m, A _1 m)} . (28) 

{m(c 3 )E%\dm=0} 

Since monopole prefactor C m is independent of A, we have immediately that 
monopoles do not contribute to the photon propagator pole. We conclude, 
in agreement with the numerical results of Section 2, that the interaction of 
the integer valued 1-form q and the gauge field A generate the photon mass. 
In other words, as it follows from Eq. (25), the Dirac strings "without" 
monopoles are responsible for the photon mass. 
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Appendix A Landau and cLandau Gauge 



Gauge-fixing is determined by extremizing with respect to gauge transfor- 
mations V x = exp{— i(f>(x)} the action [5, 6] 

s L = £ s^ = -J2 cos(Mx) - <f,( x + /*) + ^x)). (A.t) 



X [1 = 1 



At a site S% is extremized exactly by choosing <f>(x) = — arctan(a/6) where 

D 

a(x) = ^ sin A^x) — sin A^x — /i), (A. 2) 

[i=i 

D 

b(x) = ^ cos A^x) + cos — fi). (A. 3) 

Arctan is multivalued. Landau gauge is given by the branch choice 

f[-M) b>0; 
Landau gauge : </> £ < [f , 7r) 6 < 0, a < 0; (A. 4) 

([-7r,-f) 6<0, a>0; 

which minimizes This branch choice requires simultaneously that a(x) = 
and — > in the /3 — s- oo limit, and agrees with the Landau gauge of 
Ref. [6]. 

The Landau gauge used in Ref. [12] is not this but an inequivalent exten- 
sion of continuum Landau gauge. To distinguish it from the Landau gauge 
in this paper and Ref. [6], we call it "cLandau" gauge in this paper. The 
branch choice for cLandau gauge is 

7T 7T 

cLandau gauge : <f) £ [— — , — ) V a, b. (A. 5) 

This branch choice also converges to a(x) = but, unlike Landau gauge, 
preserves the discrete global "compactness" symmetry 

A " ^^-(tt + AA ifA„<0. {AM} 
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In the f3 — > oo limit ol cQED in cLandau gauge, hall ol the tends to 0, a 
quarter tends to tt and the remaining quarter to — tt. 
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